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ABSTRACT 


2 2 

The motions of thin and thick low density (q = a>p /a> c « 1) crossed- 
field electron beams have been studied on the computer. Theoretical indi- 
cations that such beams can be stable or unstable depending on the geo- 
metrical and electrical arrangements are supported, and the large amplitude 
behavior of a typical unstable configuration is found. In particular, thick, 
re-entrant beams can be quite stable against perturbations of the diocotron 
type. 
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I. INTRODUCTION 


Crossed-field electron beams (i. e. , electrons traversing crossed 
electric and magnetic fields) are characterized by the non-dimensional 
ratio 


q = oo ^/u) ^ - nm/€ n B^ (1.1) 

It has been shown'*’ that for thick, re-entrant beams having low values of 

q( < l/20), a substantial degree of stability is possible; that this stability 

can be achieved in practice is strongly implied by empirical considerations 

2 

related to the success of the Vac-Ion Pump. Two new applications for low 

3 4 5 

density crossed-field electron beams, * * have stimulated interest in this 

1 6 7 8 9 

stability, and have given rise to new theoretical studies, 9 9 9 9 experi- 

ments, ^ and to the computer experiments which are the subject of 
this paper. 

For the calculations reported in this paper, we have chosen initial 
conditions susceptible of analysis, at least in the linear regime. The agree- 
ment between the computer calculation and the analysis in the linear regime 
gives us confidence that more complicated initial conditions can be chosen, 
and their development in time followed on the computer. In addition, our 
computer results in the non-linear regime are quite informative. 

The results of the computer calculations support the view that thick, 
re-entrant low density crossed-field beams are stable; in a case for which 
the linear theory predicts instability, growth is indeed observed, but the 
mode saturates without loss of electrons to the walls. 
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II. THE MODEL 

The model used in our study is exhibited in Fig. 1. A beam of 
electrons of uniform number density n and width 2a is located sym- 
metrically between two conducting plates separated by a distance 2b(> 2a). 
A uniform magnetic field of strength B points in the z-direction. The 
model is two-dimensional, no perturbations in the z-direction being con- 
sidered. In view of this assumption, any potential applied between the 
conductors can be removed by transferring to a frame moving in the x- 
direction with a suitable velocity. Such a transformation can have no 
effect on the presence or absence of stability, or on any motion of the 
electrons in the y-direction. We can therefore consider both conducting 
plates to be at ground potential; the potential depression caused by the 
electrons in the unperturbed state is given by: 

2 a (b - y) a < y < b 

2ab - - y^ - a < y ^ a (2. 1) 

2a(b + y) - b < y < - a 

This potential depression represents the mutual repulsion of the electrons; 
the basic tendency of the electron beam is to spread out in the y-direction, 
with the electrons arriving eventually at the conducting boundaries. 

In the steady state, the divergence of the electric field in the beam 
gives rise to a shear in the drift velocity. The magnitude of this shear is 
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Fig. 1 Schematic diagram of the model chosen for the computer calculation. 
The model is two-dimensional (no variations in the z-direction) . 

It is also, as shown, periodic in the x- direction with maximum 
wavelength equal to the distance in the y- direction between the con- 
ducting electrode boundaries. The computer calculation is thus 
restricted to a square of side Zb. In the frame shown (center of 
the electron beam at rest), both electrodes are grounded, and the 
potential represents the mutual repulsion of the electrons. Note the 
shear in the electron drift velocity. 
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ne 


( 2 . 2 ) 




^< e / b > = *$■ 


= u) / a) = qw 
p ' c c 


The shear is therefore independent of the electron mass, and much smaller 
than the gyro frequency. It is assumed that the quasi-static approximation 
is valid, that is, that E/B «c . 

We can form a characteristic gyro-radius by considering the elec- 
tric field at the edge of the beam E = nea/eQ . The gyro-radius at the 
drift speed appropriate to this electric field is: 


m(E/B) 

eB 


mn 


<0 B 


qa . 


( 2 . 3 ) 


Thus for q « 1 , and for perpendicular electron thermal velocities Vj^ on 
the order of E/B , the characteristic gyro-radius is much smaller than 
the beam width. 

As stated in the Introduction, this model is susceptible of analysis 
in the linear regime. Thus, before discussing the computer calculations, 
we briefly describe the linear analysis, thereby indicating what is to be 
expected from the computer. 
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III. THEORETICAL CONSIDERATIONS 


The theoretical treatment of the small amplitude motions of 
crossed-field electron beams that follows makes use of an Eulerian treat- 
ment. This treatment is justified when, as in this case, we are studying 
wavelengths much greater than the characteristic gyro radius. 

In the limit of small q (small electron mass, low density, high 
magnetic field) the equations of motion of the electrons can be approximated 
by the equation 


E + v x B = - V(j) + v x B = 0 , (3. 1) 

This is essentially the guiding center approximation. Since the three 
vectors in (3. 1) are mutually perpendicular, 


v = - V<f> x B/B 2 . (3. 2) 

rw / “ s- ' 

Notice that the velocity vector of equation (3.2) is solenoidal. The diver- 
gence of equation (3. 1) yields 


n 


e o B 


[ curl V ] 



The equation of continuity 


Dn 

DF 


e o B 


e 


D 

Dt t curlv l 


0 


(3. 3) 


(3.4) 
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completes the system of equations. At the conducting boundaries 

= V B = 0 (3. 5) 

Since v is solenoidal, it can be considered to represent the flow of an in- 
compressible inviscid fluid. Similarly, Eq. (3. 3) shows that the electron 
density is analogous to the vorticity. Eq. (3.4) is then the conservation of 
vorticity, and Eq. (3. 5) the condition of no normal flow at the boundaries. 
Hence a complete analogy exists between the low density-crossed field elec- 
tron beam, and the sheared flow of an incompressible inviscid fluid. The 
stability of such beams is then analogous to the Kelvin -Helmholtz instability 
of sheared fluids. ^ In the case of electron beams the instability is 

usually called "slipping stream" or "diocotron. " For perturbations pro- 
pagating in the x -direction according to f(y) exp|i [kx - cot]} (k real), the 
only waves of interest are two surface waves, associated with the two free 
edges of the beam. For these two waves one obtains the dispersion relation: 

fe) 2 ■ [^-( ch tv.a 2ka )] 2 - [4^] 2 < 3 - 6 > 

This dispersion relation is shown in Fig. 2, for b/a = 1, 2, 3 and 8. 
The principal features of the dispersion relation are as follows: when the 
right hand side of Eq. (3. 6) is positive, the two waves propagate in opposite 
directions with real phase velocities of equal magnitude. When the right- 
hand side of Eq. (3. 6) is negative, the two waves have equal and opposite 
pure imaginary phase velocities, so that one wave represents pure exponential 
growth and the other pure exponential damping. This symmetry comes 
about as a result of the symmetric nature of the model in the chosen 
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I m (o)/(n 0 ) 
Re (c o/w 0 ) 


b/a = 2j 


r b/a = 8 


FUNDAMENTAL 
FOR b/a = 8 


SECOND HARMONIC 
FOR b/a = 8 


The dispersion relation (Eq. 3. 6) appropriate to the configuration 
of Fig. 1. For b/a < 2, stability is always present, but for 
b/a > 2 sufficiently long waves are always unstable. The curves 
for b/a = °o are indistinguishable from those shown for the case 
b/a = 8. When the wavelength is restricted to take on the values 
2b, or sub-multiples of 2b, the abscissa is "quantized. n The first 
two permissible modes for the case b/a = 8 are shown; only the 
fundamental is unstable. 
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frame of reference. For b/a < 2 , stability is always present, while 
for b/a > 2 , long waves are unstable, up to some critical wave number; 
beyond this wave number, stability is present. 

It was stated in the Introduction that a condition for stability of 
this type of crossed-field electron beam is that the beam should be re-entrant. 
The reason for this can be seen from Fig. 2, for if arbitrarily long waves 
are permitted, beams for which b/a > 2 will always be unstable. In 
practice, the re-entrant character of the beam is obtained by using cylin- 
drical or toroidal arrangements. For our present model, the re-entrant 
feature is simulated by requiring the beam to have a finite periodicity in the 
x-direction. In particular, we require all potentials to be periodic in the 
x-direction with periods equal to 2b and sub-multiples of 2b . Thus the 
computer calculation is restricted to a square area of side 2b . Three 
squares of this kind are shown in Fig. 1. From the point of view of Eq. (3. 6), 
we require k to have only the values 

ka = m 7r , m = 1 , 2 , . . . (3.7) 

For the case b = 8a , these wave numbers are marked on Fig. 2, for 
m = 1 and 2; it will be seen that only the fundamental mode m = 1 is long 
enough to be unstable, with a growth rate equal to 0. 2 Wq . 

The inertial terms ignored in Eq. (3. 1) are characteristically 
smaller than the electric and magnetic terms by the factor q . The 
inclusion of these terms has been shown (in a cold approximation) to yield 
a short-wave instability with a growth rate ~ co^e . In practical cases, 
this growth is negligible for q ~ l/20 . 
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The linear theoretical analysis gives no information on the behavior 
of the system after saturation. However, the aerodynamic analogy suggests 
the formation of a "vortex street. " The formation of such a street is shown 
by the computer calculations (Section 5). 
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IV. COMPUTATION TECHNIQUE 


In the computer calculation, the continuous distribution of charge 
implied by Eq. (2. 1) is replaced by 1,000 discrete "rods" of charge whose 
motion is calculated step-wise in time subject to a self-consistent electro- 
static field. This is essentially a Lagrangian technique. Such a technique 
is, of course, perfectly sound; however, unlike many neutral plasma and 
high density (q ~ 1) electron beam calculations, it is not positively 
required. This point was made in justification of the Eulerian approach to 
the analysis of the previous section. 

For our calculation, the Lagrangian approach was chosen because 

1 5 

of the existence of a substantial program in this form, developed by one 
of us (R.W.H.) for use in the study of the anomalous diffusion of plasmas. 
This program required only slight modification before it could be used for 
the present studies. 

The positions and velocities of the 1,000 charged rods are advanced 
stepwise in time as follows: (1) from the position of each charge, one 

derives a charge distribution; (2) from the charge distribution, using 
Poisson's equation, one calculates the electrostatic potential; (3) from the 
potential, one derives the electric field acting on each charge; (4) one then 
integrates Newton's law of motion for short time step DT, using the local 
electric field and a constant external magnetic field. This calculation gives 
a new position and velocity for each rod, and the cycle repeats at step (1). 
As an initial condition the 1,000 rods are distributed randomly within a band 
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in the interior. The initial velocities are selected to give a Maxwellian 
distribution with the temperature as an input quantity. The errors inherent 
in such a model are discussed in Ref. 15, and the special rapid technique 
used for the solution of poisson's equation is described in Ref. 16. 

An interesting feature of the present calculation arises in connection 
with step (4), the integration of Newton's law of motion. The full equation 
of motion 


d^x T dx ”1 

m — j ~ ~ e E + — x § > (4. 1) 

dt L dt J 


is represented by the central difference approximation 


[ x(t + DT) - 2x(t) + x(t - DT)] 

m > 5 

DT 


e E(t) + 

* 


jx(t+ DT) - x(t-DT)| 


2DT 


B 


(4.2) 


where x is the position vector of a gyrating rod. This scheme was 

1 7 

recommended by Buneman, and has the important property of being rever- 
sible in time; it can also be used for time steps greater than a gyro-period, 
provided that the net displacement in each time step is small compared 
with the scale size for changes in the electrostatic potential. To see this 
we note that the solution in constant E and B fields to both the differential 
Eq. (4. 1) and the difference Eq. (4.2) can be expressed as a gyration with 
Larmor radius R and gyro -frequency <x> c about a guiding center drifting 
with a velocity v^ . That is to say, taking axes such that E x = 0 , 
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(4. 3) 


x(t) = Vj-jt + R cos(oi> c t + <j) ) + constant 
y(t) = R sin (w^t + $) + constant 


where <j) is a constant phase determined by the initial conditions. 

It is easy to relate the parameters of the motion derived from the 
difference equation (distinguished by a star) to the parameters of the exact 
motion derived from the differential equation. For simplicity we take the 
case that the electron has initially zero velocity (the cold case of Section V) 
and obtain the relations: 



a- •y 

go = arctan (oo^DT/2) • 

R* is the radial distance of the vertices of the polygon which the electron 
follows in this numerical approximation. It differs as indicated from the 
true gyro -radius R . is the frequency with which this polygon is 

traversed. This frequency differs as indicated from the true gyro -fre- 
quency . We note that the difference equation reproduces exactly the 
guiding center drift. The only errors introduced by differencing are in the 
gyro -frequency and Larmor radius. 

For the study of the linear growth regime in Fig. 3, covering five 
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LOGARITHM OF THE AMPLITUDE 
OF THE FUNDAMENTAL MODE 



Fig. 3 Growth in amplitude of the unstable fundamental mode for the thin 

beam (b/a = 8) case. The initial amplitude results from the random 
initial position of the electrons in the beam, and the growth rate is 
derived from Fig. 2. The time step used for this calculation was 
about l/5 of a gyro period, or l/lOO of the shear period Ta = 27 t / u ) q . 
The indicated exponential growth rate is approximately 0. ZcOq. 
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growth times, we have been able to use five time steps in a cyclotron 
period. This leads to 

R r = 1.1 8R and to = 0.89 to 

c c 

The angle turned through between samplings of the electric field is 
co^DT = 66°. With these values one has a reasonable representation of 
the full dynamic motion of the electrons. 

In the study of the long term non-linear behavior in Figs. 4 and 5. 
covering 25 growth times it was necessary to use a longer step which was 
approximately equal to the gyro -period. In this case 

R* = 3. 3R = 0. 52 V d DT 

a* 

go * = 0. 4 go and 

c c 

5jc O 

go DT = 144° . 
c 

The Larmor radius and gyro period are now badly distorted, hence such a 
large step would not be acceptable in the study of effects which depended on 
electron inertia. However, since we are studying the diocotron instability 
which is independent of the electron mass, this is unimportant. It is only 
necessary that the correct guiding center drift motion be obtained. This 
will be the case if the distance moved by the guiding center in a time step 
and the Larmor radius are small compared with the wavelength, X , of the 
diocotron wave. The relations (4.4) show that for large time steps the first 
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'ig. 4 Frames from the movie exhibiting the growth and development of the 
thin unstable electron beam. The time step used in this calculation 
was approximately one gyro period, or l/ZO the shear period Tq. 
Times on this figure are in units of 0. 98 tq. The relatively long 
time step used in the calculation should not affect the long wave 
behavior exhibited in this figure. 
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E-FIELD ENERGY (ARBITRARY UNITS) 



Fig. 5 Energy spectrum associated with the computer run exhibited in 

Fig. 4. The unstable fundamental mode grows and then saturates. 
The total energy remains constant to within 3-l/2%, The energy 
residing in the higher modes starts out at about the same level as the 
fundamental mode (about 10“^) and does not increase significantly 
during the run. 
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condition guarantees that the second is satisfied. In our case we have 
taken q = which leads to 

V n DT R* 

— ^ — — 0 . 02 , = 0 . 01 , 

where X = 2b is the wavelength of principal interest. The description of 
very short waves is, however, likely to be in error. The shortest wave- 
length available to the program is ~ 2b/24 . The sampling of the field 
every 144° ensures that the field is evenly sampled on either side of the 
line of motion of the guiding center. Thus we expect the guiding center 
motion to be correctly represented even though the time step is too large 
to represent the inertial aspects of the electron motion correctly. 

The main output from the computer model is produced on a Strom - 
berg-Carlson 4020 computer recorder. This maybe a complete record of 
the rod positions and potential distribution or alternatively simply a record 
of the rod positions. The latter form of output is much less costly in 
computer time and is shown in Fig. 4. Either of these displays can be 
brought out on movie film for projection purposes. A time step, with the 
record of rod positions only, takes 2. 1 secs on an IBM 7090 computer; the 
total computer time for a typical run of 250 steps with movie output was 
17 minutes. 
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V. RESULTS 


Linear Growth Regime 

Two cases have been run, namely b/a = 2 (the thick stable beam) 

and b/a = 8 (the thin unstable beam). For both cases q was taken as 

1 8 

l/20, which is the same value used by Byers in a similar experiment. 

This is sufficiently small to ensure that we are studying the long wavelength 

^ j 

phenomena, since the short wave growth rate is estimated to be c^e ~ 
18 

4x10 gOq - leading to a growth time roughly 10 times longer than the 

calculation time. On the other hand q is sufficiently large for us to take 

approximately five steps of the calculation in a gyro period of the electron. 

Five steps per gyroperiod has been used successfully previously by 
1 5 

Hockney and enables the full orbit of the gyrating electron to be followed. 
These results thus include the inertia of the electrons and we do not make 
the approximation of guiding center motion. 

For the long wave diocotron instability the growth time is of the 
order of Tq = 27t/g0q (see Fig. 2). In the case of the thick beam there was 
no measurable growth of any instabilities over a period of five growth times, 
and the beam remained well confined to its initial width. 

In the case of the thin beam we find, as predicted theoretically, that 
the m = 1 wave is unstable. Fig. 3 shows the excellent agreement that 
was obtained between the theoretical growth rate and that observed in the 
computer model. This agreement gives us confidence in believing the new 
results which are found with computer model and which are not predicted 
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by linear theory. These relate to the large amplitude, non-linear behavior 

of the system. We observe in Fig. 3 that the amplitude of the m = 1 wave 

saturates after about 3 Tq . At this time the beam has twisted itself into a 

vortex but the amplitude of the vortex is not sufficient to drive any of the 

1 8 

electrons to the wall. This result is in agreement with that of Byers who 
found that the Kelvin -Helmholtz instability gave rise to confined stable 
vortices in a two-fluid Eulerian model of a plasma. 

The above runs were performed for two initial conditions. In the 
first, or cold case, the electrons, which are given random positions within 
the beam, have zero initial velocity. The kinetic energy which they quickly 
acquire is due entirely to the electrostatic field of the beam itself. In the 
second case a small initial velocity, selected from a Maxwellian distribution, 
is given to the electrons. The characteristic thermal velocity V ^ = 

i /z 

(kT / m ) v ^ was taken as approximately one half of the maximum shear 
velocity of the beam. We observed no difference in the macroscopic long 
wave behavior of either the thick or thin beams. This result is in agree- 
ment with the fact that the long wave instability is predicted from guiding 
center theory which is not influenced by the electron thermal velocity. 

Long Term Non-linear Behavior 

For the long term behavior we have used a time step approximately 
equal to the gyroperiod, and the results have been displayed in a short 
16 mm movie. Still frames from this movie are reproduced in Fig. 4 in 
which each electron-rod is represented by a dot. As expected from the 
analogous situation in hydrodynamics, the beam breaks up into a "vortex" 
street. The vortices are completely formed at t = 5. OTq and persist 
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without change for at least the next twenty shear periods. 

As a check on this long term behavior we have also computed the 
variation of electric -field energy during this time, and its division between 
the different wave numbers. This is shown in Fig. 5. We find that the 
total field energy drops by 3-l/2% during the 500 time steps represented 
in Figs. 4 and 5. This represents the cumulative effects of all the finite 
difference errors arising in the computer model. These are the finite 
size of the space -mesh used in the solution of Poisson's equation, and the 
finite size of the time step. The energy in the m = 1 unstable wave is seen 
to increase by 3 orders of magnitude during the first three growth times 
and then to saturate with about one eighth of the total field-energy. The 
energy which goes into the m = 1 wave is taken from the electrostatic 
energy of the initial beam which is the energy in the m = 0 "wave, 11 An 
interesting feature of the computer results is that this interchange of energy 
between the m = 0 and m = 1 waves persists on a smaller scale for the 
duration of the run. One can see a periodic exchange of energy with a 
period approximately equal to 2 Tq . The amplitude of this oscillation 
appears to decrease by a factor two during the run. The effect of the 
oscillation on the form of the vortex is minor and can be seen by comparing 
the vortex in Fig. 4 and t = 5 Tq when the m = 1 wave is a maximum and 
with that at t = IOTq w ^ en the wave i s a minimum. 
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VI. CONCLUSIONS 


A Lagrangian computer model consisting of 1,000 electron rods has 
been used to investigate the long wave diocotron instability in a crossed- 
field electron beam. A thick and thin beam were found to be stable and 
unstable respectively, as predicted by theory. Excellent agreement was 
obtained between the theoretical growth rate and that observed in the model. 
The computer model was able to follow the non-linear behavior of the 
system. The unstable beam broke up into a vortex street which remained 
for a long time. The size of the vortex was not enough to transport elec- 
trons to the wall of the system and the beam remained confined to the 
interior. The results of this paper have been made into a 4-minute, 16mm 
movie. 
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